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Abstract 

We derive and discuss the equations of motion for the condensate and its 



> 

co , 

<^> fluctuations for a dilute, weakly interacting Bose gas in an external potential 

within the self-consistent Hartree-Fock-Bogoliubov (HFB) approximation. 
Account is taken of the depletion of the condensate and the anomalous Bose 
correlations, which are important at finite temperatures. We give a critical 
analysis of the self-consistent HFB approximation in terms of the Hohenberg- 

o 
o 



Martin classification of approximations (conserving vs gapless) and point out 
that the Popov approximation to the full HFB gives a gapless single-particle 

x; 

spectrum at all temperatures. The Beliaev second-order approximation is 
discussed as the spectrum generated by functional differentiation of the HFB 
single-particle Green's function. We emphasize that the problem of deter- 
mining the excitation spectrum of a Bose-condensed gas (homogeneous or 
inhomogeneous) is difficult because of the need to satisfy several different 
constraints. 

PACS numbers: 03.75.Fi, 05.30.Jp, 32.80Pj, 67.90,+z 
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I. INTRODUCTION 



The recent observation of Bose condensation in atomic gases trapped in an external 
potential well has focussed attention on the excitation spectrum of an inhomogeneous weakly 
interacting Bose-condensed gas at finite temperatures. Until very recently, this excitation 
spectrum was only addressed indirectly in the course of calculating P,R| the local density 
n(r) and local condensate density n c (r). 

The excitation spectrum of a Bose system is very dependent on the subtle dynami- 
cal correlations induced by the Bose condensate M. In particular, one must calculate the 
single-particle and density-fluctuation spectrum in a consistent manner to ensure that they 
are identical, as they must be in a Bose-condensed system. In general, self-energy en- 
ergy approximations which lead to a single-particle spectrum without an energy gap in the 
long-wavelength limit are not consistent with a density-fluctuation spectrum which satisfies 
conservation laws (or the corresponding sum rules). 

A major purpose of this paper is to point out that a useful way of understanding various 
approximations for the excitations in a spatially non-uniform Bose gas is provided by the kind 
of analysis developed in the early 1960's for a spatially uniform Bose gas and codified in the 
well-known paper by Hohenberg and Martin [3J]. It is based on distinguishing "conserving" 
vs "gapless" approximations and the key role played by the Hugenholtz-Pines theorem. 

We use this classification to examine the equations of motion for a non-uniform con- 
densate and its excited states in the finite temperature Hartree-Fock-Bogoliubov (HFB) 
approximation, as well as simpler approximations discussed in the literature. An approxi- 
mate version of the HFB due to Popov |5||| gives a simple gapless approximation for the 
single-particle spectrum below the transition (we do not discuss the critical region very near 
Tbe) but also reduces to the T = Bogoliubov approximation discussed by Fetter [0] for 
a non-uniform gas. We also sketch the procedure |||| for using the HFB approximation 
for the single-particle self-energies to generate by functional differentiation a density fluc- 
tuation spectrum which is gapless in the long wavelength limit. This can be shown to be 
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identical to Beliaev's second order approximation for the single-particle spectrum || at all 
temperatures. 

The emphasis in this paper is on the formal structure of various kinds of approximations 
and how to assess their validity, using the HFB approximation as a specific example. Our 
results in (|8p and ( P3"D should provide a platform for future numerical calculations of the 
HFB excitation spectrum of an inhomogeneous Bose gas at T ^ 0, of the kind recently 



discussed by Edwards et al. [|lt|] an d Fetter [fLl[] using the Bogoliubov approximation at 



T = 0. However, as we emphasize in the final section, the question of what are the 
"correct" excitations in a dilute weakly-interacting Bose gas is not simple one. It requires 
the full apparatus of many-body theory. 

II. HARTREE-FOCK-BOGOLIUBOV APPROXIMATION 



Our starting point is the exact Heisenberg equation of motion |T2| , [T3|j for the Bose field 
operator -0(r) 

.dip(r,t) ( V 2 x 



dt \ 2m 



+ U ex {r) - fi t) + gif)\r, t)^{r, t)^(r, t) , (1) 



where we have assumed a short-range interaction v(r — r') = g5(r — r') between the atoms 
(In the s-wave approximation, which is adequate for the very dilute gases of interest, one 
has g = 47ra/m). We assume the atoms are trapped in a static external potential U ex (r). 
Separating out the condensate part in the usual fashion we have 

^(r,t) = $(r)+^M) , (2) 

where $(r) = (ip(r,t)) = plays the role of a spatially-varying macroscopic Bose field. 

The thermal average in (JJ) and elsewhere is the usual anomalous average appropriate to a 
Bose broken symmetry Using (0), the interaction term in (|I|) can be written out in the 
form 

_ |$|2 $ + 2|$| 2 ^ + $ 2 ^t + + 2$^^ + ft<ijnjt . (3) 



We treat the last term in (|3|) in the self-consistent mean-field approximation, namely 



ft{r,t)4>{r,t)iP(r,t) ~ 2(ft{r)4>{r))i>{r,t) + (^(r)^(r))^(r, t) 



(4) 



and then (|3[) reduces to 



^ t (r,t)^(r,t)^(r,t) = |$(r)| 2 $(r) 

+ 2 [|$(r)| 2 + (^t( r )^( r ))] ^( r ,t) + [$ 2 ( r ) + (^(r)^(r)] ft(r,t) 

+ 2$(r)^ t (r, t)^(r, t) + $*(r)^(r, t)^(r, t) . (5) 

For the zero-range interaction we are considering, the Hartree and Fock (exchange) terms 
are identical. This is the origin of the factor of 2 in (|5|) and subsequent equations. 

The time independent, spatially-inhomogeneous Bose order parameter $(r) is given di- 
rectly by taking an average over ([!]), 



- — + U ex (r) - n J $(r) + g(ft{T)$(T)j>(T)) = . 



(6) 



The time- dependence cancels out in the last term using the cyclic invariance of the trace. 
Taking the anomalous average of ([5]), the linear terms in ip(r,t) vanish since (ip(r,t)) = 
= and we are left with 



(ft(r)4)(r)i)(r)) = |$(r)| 2 $(r) + 2$(r)(^ t (r)^(r)) + $*(r)(^(r)^(r)) 



(7) 



Using this in @, we find 



' V 2 

-— + U ex [ 
, 2m 



r) - /i $(r) + g[n c (r) + 2n(r)]$(r) + #m(r)<I>*(r) = 



where we have introduced the local densities: 



(8) 



n c (r) = |$(r)| 2 
n(r) = (^(r)^!-)) 
m(r) = (ip(r)ij)(r)) . 



(9) 
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We note that (|8|) only reduces to a closed (non-linear) equation for $(r) when both n(r) and 
m(r) are neglected (i.e. in the Bogoliubov approximation 0). In the case of a spatially- 
uniform Bose gas (U ex (r) — > 0), all the functions in (||) become constant, independent of 
position. 

The excitations of the condensate are described by ip(r,t) in (g) and are given by the 
exact equation of motion 



.dijj(r, t) 

' dt 



2m 



+ U ex (r) - fi 



+ g\ft(r,t)j>(T,t)$(T,t) - ($ + (r)$(r)$(r)) 



. 



(10) 



This follows from subtracting (|6]) from (|l|). Consistent with our mean- field derivation of 
the quadratic terms in the last line of fl5|) are treated in a mean- field approximation , i.e. 
we use 



(11) 



With fllTD, the last term in ( pT0| ) is given by 



^(r,t)^(r,t)^(r,t) - (^(r)^(r)^(r)) 
~ 2(^ t (r)^(r))^(r,t) + (^(r)^(r))^(r, t) 



(12) 



Using (1121) , ([TOP reduces to 



2m 



+ [/ ex (r)-// U(r,t) 



+ 2gn(r)ip(r, t) + gm(r)ip'(r, t) , 

where the self-consistent densities are defined by (see (0)) 

n(r) = (^(r)^(r)) = |$(r)| 2 + (^(r)^(r)) 
= n c (r) + n(r) 



(13) 



(14) 
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m(r) = (if)(r)if)(r)) = $ 2 (r) + (ip(r)ip(r)) 

= $ 2 (r) +m(r) . (15) 

One can easily derive the analogous equation of motion for if)'(r,t). 

The coupled equations of motions given by (H) and (p~3| ) correspond to the Hartree-Fock- 
Bogoliubov (HFB) approximation. One can derive the same or similar results by a variety 
of ways but the above approach has the advantage when one is interested in the excitations 
since it is easy to relate these equations of motion to a more general Green's functions [§,§3 
formulation discussed in Section III. 

The preceding derivation of (|8|) and ([13]) is equivalent to the grand canonical Hamiltonian 

K = H-/iN = J drft(r) 

+ - J dr^ + (r)^ + (r)^(r)^(r) (16) 

being treated using a self-consistent quadratic approximation. Expanding the field operators 
in this expression using (H), one finds after a little algebra |L4j] 

Khfb = J dr$*(r)[f(r) + ^n c (r)]$(r) 
+ J drft(r)£iP(r) 
+~ J drm(r)if)* (r)if) + (r) 

+~ J drm*(r)^(r)^(r) , (17) 

where £ = T(r) + 2gn(r), with 

f{r) = ~ + U a (r)-ii. (18) 

The coefficients of the linear terms in ip, ip* in ( |17|) can be shown to vanish by using the fact 
that 3>(r) is given by the solution of @. In going from (jlBf ) to (|17|) , the on/?/ approximation 
involves how we treat the terms cubic and quartic in if) and if)*, namely (all quantities depend 
on r) 



V 2 
2m 



fi 



1>{t) 
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K 3 = g I dr^ft^ + $ ft flip] 

~g J dv[2n<$>* + m*$]^ + g J dv{2n§ + m§*]ft (19) 



= — I dvftftipip 
9 

~ 2 



dv[4nftip + m*ip-ip + mftft) . (20) 



That is, within the HFB mean-field approximation, K% is linear in ip, ft while is quadratic 
in ip,ft. We note that the coefficients in ( |T9"D and (|20| ) involve either n(r) or m(r). 
One can easily diagonalize (|17D by using the linear transformation 

j 

ft(r) = Y.{n){v)a]-v ] {r)a ] ) , (21) 

i 

where a-j and a] are annihilation and creation operators satisfying the usual Bose commu- 
tation relations. One can show that (|17]) reduces to ]15| 

K HFB = J dr$*(r)[ f (r) + ^|$(r)| 2 ]$(r) 

-E S i/*l v i( r )l 2 + E^«]«i ( 22 ) 
i J j 

if the c-number functions Uj(r) and fj(r) are given by the solutions of the coupled eigenvalue 
equations 

tuj{v) - gm{r) Vj {r) = E jUj {r) 

Cvj{r) — gm*(r)uj(r) = —EjVj(r) . (23) 

Here the operator L is defined before ( |TBD and m(r) is defined in (|T3|) . One can also derive the 
same results working directly with the equation of motion (O) and its Hermitian conjugate. 
In particular, one may easily verify that 



-iEjt „,* /_\ * \AEjfr\ 

^M) = E[«J(#^ - "iW^ ( 24 ) 



solves (|I3D (and its adjoint) if Uj,Vj and £7j satisfy the generalized "Bogoliubov" equations 
given in (^). The results in (|22| ) and (^) effectively reduce the problem to a gas of non- 
interacting Bose quasiparticles with an energy spectrum given by Ej. 

We note that with (pTD, one easily can derive expressions for n(r) and rfi(r) in (|9|) in 
terms of the self-consistent solutions of the coupled equations (|2"B|), namely 



n(r) = E { [k«| 2 + KMI 2 ] WO + I^W| 2 } 

3 

m(r) = -E«,(rK(r)[2iVo(^) + 1] , (25) 

3 

where N (E) is the Bose distribution for the quasiparticle excitations 

<$ & i) B FB = = N o(E,) ■ (26) 

One must solve the coupled HFB equations (|23|) and the condensate equation (§) using self- 
consistent values of n(r) and m(r), given by (|25|) , and n c (r) = |$(r)| 2 . Computationally, 
solving this set of coupled equations is similar to solving the simpler T = Bogoliubov 
approximation The major difference is that (H) is no longer a closed equation 

for $(r) but involves the self-consistent diagonal and off-diagonal densities, n(r) and m(r). 
The latter depend on the excited states given by the solutions of the self-consistent coupled 
equations in (f23|). 

It is useful to summarize how one would solve the HFB equations given by (§) and ( p5| ) . 
for a given anisotropic parabolic potential-well trap and temperature: 

(a) First calculate rio(r) = |$(r)| 2 and n(r) for a non-interacting gas (g = 0). In this case, 
the off-diagonal correlation function m(r) = 0. 

(b) Use these non-interacting gas results in ( p3|) to solve for Uj, Vj, and Ej. With these 
excited states, one can calculate n(r) and m(r) from (^) and use these in (§) to find 
$(r). 

(c) Iterate the above procedure to self-consistency. 



S 



Before discussing the implications and validity of (§) and fll3|), it is useful to relate these 
results to simpler approximations given in the literature on inhomogeneous dilute Bose gases. 
We recall that ( |T3| ) and its adjoint are equivalent to (^3|) and fl24|). If we set both n(r) and 
m(r) to zero, (|8]) reduces to the well-known Gross-Pitaevskii approximation, where <J>(r) is 
given by a closed non-linear Schrodinger equation (NLSE) |16|| . Neglecting n(r) and m(r) 
in (|l~3| ) or ( p3|) leads to the standard Bogoliubov approximation. This simplified version of 
Eqs. (H) and (|13|) has been derived and discussed in detail by Fetter |7j (see also Refs. 
and Since one is ignoring the non-condensate atoms (n = 0), this approximation is 

only appropriate at T ~ 0, where most of the atoms are in the condensate (n c ~ n). Note 
that in the Bogoliubov approximation, the cubic terms (Ks) and quartic terms {K4) in ([19]) 
and fl20|) are omitted completely. 

The next level of approximating (|8|) and (^) would be to keep n(r) but neglect the 
anomalous density m(r) as being small compared to both n(r) and no(r). This approxima- 
tion has been used by Popov || (in a homogeneous gas) to discuss the finite-temperature 
region close to the BE transition. This Popov approximation || formally reduces to the 
Bogoliubov approximation at T = 0, where n also becomes negligible. As we shall see in 
Section III, it leads to a gapless spectrum, in contrast with the full HFB. Thus omitting m(r) 
but calculating n(r) in a self-consistent way seems to give a reasonable first approximation 
for the excitation spectrum in Bose gases at all temperatures. 

The calculation of n(r) and n c (r) given by Goldman et al. involved equations equiva- 
lent to (§) with m(r) = and to (^) with m(r) = 0. This somewhat ad-hoc approximation 
is equivalent to putting Vj(r) = in (|23|). This simplified version of (p3|) has been justified in 
Refs. [0,|J as being adequate for the purpose of calculating n(r) and no(r) at temperatures 
just below Tbe-i even if it does not lead to a very good approximation for the excitation 
spectrum. 
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III. CONSERVING VS GAPLESS APPROXIMATIONS 



We now turn to a discussion of the HFB equations given by @ and (0). The general 
problem of finding the excitation frequencies in a homogeneous weakly interacting Bose- 
condensed gas was exhaustively studied and various approximations were classified in a 
famous paper by Hohenberg and Martin (see in particular Section VI of Ref. 1J). Any 
approximation put forward in the study of an inhomogeneous Bose gas can be usefully 
analyzed by examining its implications in the corresponding homogeneous case. 

For interacting Bose-condensed systems, it is convenient to formulate the discussion in 
terms of Green's functions [@,@. A key role is played by the 2x2 matrix single-particle 
Green's function, defined as 

0,(1, 1') = -i<Tfc(l)fct(i')) (27) 

where 

/ ^(1) \ 

*(1)= ; *t(i) = (Vi t(i) )Vi (i)) (28) 

and 1 represents r, t. Writing out the 2x2 matrix in (|27|) explicitly, we have 



/ i(T^(l)^(V)) z(T^(1M1')> \ 

^(1,1') = - . + . + _ , N ■ (29) 

\i(T^ (1)^(1')) i(T^{l)^{l'))j 

Separating out the Bose condensate part of the field operator as in (|]), (|2^) naturally splits 

into two parts 

G aP (l, 1') = G aP {\, 1') + G hoi {l)G\ p {l') , (30) 
where the condensate Green's function is described by 

Gi(l) = V^<*(l)) = V=i | (31) 

\$*(1) 

61(1) = >/=<(*t(l)> = >/=i(<&*(l),*(l)) • (32) 
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and G a p(l, 1') is identical to ( p9|) except that it involves the non-condensate part of the 
field operators (if) and if)*). The signature of a Bose-condensed interacting system is the 
appearance of anomalous correlation functions such as $(1) = (V>(1)) and 612(1, 1') = 
-i(T^(l)^(V)). 

A very convenient way of generating the equations of motion for G\(l, V) and Gi(l) is 
to use functional differentiation with respect to external fields Q. For homogeneous systems 
in thermal equilibrium, it is customary to set these generating fields to zero at the end of 
the calculation. However, in the context of the newly observed Bose-condensed gases [|TJ, 
these external fields are of direct physical relevance, playing the role of the static external 



trapping potential as well as time-dependent external perturbations |10| J11 |. In this regard, 
the general formalism developed in Ref. Q is particularly useful. The general equations of 
motion for Gi(l, 1') and Gi(l) are given by (2.6) and (2.7) of Ref. ||. If we only include a 
static external potential, the equation of motion for G\(l, 1') is 

[iT®^-f(v)8 a p}Ge a/ (l,l') 
= 5(1 - 1')<W + X Q/3 (1,2)GV(2, 1') , (33) 

/ 1 \ 

where = and the single particle operator T(r) is defined in (|18D. Here, 

\0 -lj 

repeated Greek indices are summed and a bar represents the usual integration over r, t. 
This Dyson-Beliaev equations in (|33|) for the single-particle Green's functions in a Bose- 
condensed system depends on the 2x2 matrix self-energy E Q!| g(l, 1') in addition to the 
chemical potential fi. The analogous equation of motion for the two components of Gi is 
given by (for a static external field) 

-f(r)G, a ( r ) = V^Va(r). (34) 

Here we have introduced the condensate "source function" rj a (r), which itself is a functional 
of G1/2 and G\. Writing (|34]) more explicitly, it corresponds to 



2 



V 

— + U ex (r) - 11 



$(r) = -^i(r) . (35) 
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One can now define what a "conserving" approximation is, namely that there exists a 
functional $[Gi,Gi] of the correlation functions Gi and G\ such that the self-energy and 
the condensate source functions are given by 



<J$[Gi,Gi] 
1 <y$[Gi,G?i] 



£(l,l') 
77(1) . 



(36) 



2v^ <JGt(l) 

2 

Martin and De Dominicis JT7J have proven that if such a functional $ exists, the two-particle 
Green's functions generated by functional differentiation of Gi(l, 1') [given by ( p9|) - (p2|)] 
with respect to a time- dependent, spatially- varying external field are guaranteed to satisfy 
the usual conservation laws [|Tj|. The HFB approximation discussed in Section II can be 
shown to be an example of such a conserving approximation. Using Eq. (6.42) of Ref. ||, 
one can easily verify that 



J HFB 



Xl') 
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2n(r) m(r) 



m*(r) 2n(r) 



5(1-1') 



(37) 



and 



n(r) + n(r) 



Vhfb(1) = 9 



(38) 



m(r) \ / $(r) 
m*(r), n(r) + n(r)J \&*{r) 
When substituted into (|33D and (p4]), one finds results equivalent to p3| ) and (||), respectively. 
This HFB approximation is often called the Girardeau- Arnowitt approximation in the Bose 
gas literature PJ. Solving ( ^3]) using (|3"7|), we arrive at the coupled equations for Gu and 
G21, 

d 



1— - T(r) - 2^(r) 
.9 

-z— - T(r) - 2^n(r) 



Gu(l,lO- W)G a i(l,l') = 5(1-1') 
G 21 (l,l , )-^*(r)G 1 i(l,l , ) = 0. 



(39) 



This is the HFB generalization of the Bogoliubov approximation, the latter being given by 
(55.22) and (55.24) of Fetter and Walecka [0 for inhomogeneous systems. As expected, the 
structure of (|39T) is identical to (E23|). 
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The problem with a conserving approximation is that while the two-particle Green's 
functions [such as the density response function Xnn(l> 1')] are guaranteed to have a spectrum 
consistent with conservation laws such as the equation of continuity ||18|| , the single-particle 



spectrum of G\(l, 1') which one starts with may have unphysical features. In particular, 
one is not ensured of a gapless spectrum for long wavelength excitations. In homogeneous 
systems, such a gapless spectrum of Gi(l,V) is guaranteed if the self-energies satisfy the 
Hugenholtz-Pines (HP) theorem, namely P,|l9| 



H = E n (q = 0, u = 0) - Si 2 (q = 0, u = 0) . (40) 

More generally, it can be proven that this theorem will be satisfied if the self-energies can 
be generated from the condensate source function according to H 

^i^M- = E(l, 1') . (41) 
SGx(V) 

2 

Here r)(l) is an explicit functional of Gi and implicitly through the dependence of G\ on 
Gi. Approximations for the self-energy consistent with ( |4"T| ) are referred to as "gapless" 
approximations. The HFB approximation for 7/(1) and E(l, 1') do not satisfy (0), and as 
a consequence the HFB spectrum of Gi(l,V) has a quasiparticle spectrum with a gap in 



is 



the limit of long wavelengths [§]]. In contrast, the simple Bogoliubov approximation [12 
a gapless approximation, as is the Beliaev second-order approximation ||. Figs. 4 and 5 
of Ref. Q give a convenient summary of conserving and gapless approximation for uniform 
Bose gases. 

We can illustrate the above somewhat formal discussion by considering the HFB approx- 
imation for a uniform system (U ex (r) — ► 0). Using (^) in (0), or equivalently (§), we have 
(since $(r) = constant) 

- ^$ + g(n + n)$ + m$* = (42) 

and hence the HFB gives the following result for the chemical potential ||, 

/j, = g{n + ft + in) . (43) 
13 



For a homogeneous system, where we have G a p(l, 1') = G a p{l — 1'), one can solve (|39| ) by 
Fourier transformation. The poles of Gapik, uj) are easily found to be given by 



2 



- A* + 2<?n) - G?m) 2 . (44) 



E, 

Making use of the HFB chemical potential \i given by fl4"5|), one sees that there is a finite 
energy gap at k = 0, 

El=o = 9 2 W ~ m] 2 - g 2 [n c + fh] 2 = Ag 2 \fh\n c , (45) 

where we recall that m is negative [see fl25f )1. Clearly the HFB does not satisfy the HP 
theorem (f40|) , since the latter relation gives 

/i = 2gn — gm = g(n + n — m) , (46) 

which differs from the HFB result in (f43D- Using (^) in (|44|) gives a gapless spectrum but 
this procedure is ad hoc. One is using a chemical potential which is not consistent with the 
equation of motion for G*i(l) or, equivalently, $(r). 

At the end of Section II, we introduced several approximate versions of the HFB results. 
We can now discuss these approximations in terms of the excitation spectrum they lead to 
in the homogeneous case: 



(a) "Bogoliubov" corresponds to setting m = n = 0. Both (H3T) and (HBT) reduce to the 



same result [i = gn c and hence the single-particle excitations are gapless. This is valid 
at zero temperature. 



(b) "Popov" corresponds to setting fh = but keeping h finite. Both (f43D and (|6|) reduce 
to /i = 2gn — gn c and hence the excitations are gapless. As discussed at the end of 
Section II, this gives a reasonable approximation at all temperatures. 

(c) In Goldman et al. the off-diagonal self-energies in ([57]) involving m(r) are neglected 



but only m(r) is neglected in the source term (|3q) . The result is that fl44|) reduces 



to Ek = — A 4 + 2gn while (|^) reduces to pi — gn + gh — 2gn — gn c . Thus one 
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obtains a spectrum Ek = + gn c with an energy gap. This approximation can 
be used for calculating n and other thermodynamic quantities near Tbe, where the 
thermally significant excitations have a large enough wavevector that the spectrum 
E k = + gn c is adequate |20 . 



In the Popov approximation to the HFB (in which the anomalous density m(r) is omit- 
ted), there is a special solution of the generalized Bogoliubov equations ( f23"D corresponding 
to Uq(t) = Vo(r) = $(r) with E = 0, where $(r) is a solution of (|8]) with m(r) = 0. As with 
the simpler Bogoliubov approximation discussed by Fetter 0,0, one sees in this case that 
the condensate wave function $(r) plays the role of the zero energy single-particle mode and 
the fluctuations of the condensate associated with ip are described by all the higher energy 
modes (Ei > 0). In contrast, within the full HFB, there is no zero energy solution of (|23|) 
corresponding to the solution $(r) of (||). 

As we have noted above, conserving approximations for G\(l, 1') and Gi(l) such as the 
HFB can be used to generate two-particle Green's functions which are guaranteed to satisfy 



conservation laws and related sum rules 0,17,18] and thus will lead to a gapless density 



fluctuation spectrum. In addition, we recall that in the presence of a Bose broken symme- 
try, the single-particle and density fluctuation correlation functions exhibit the same poles 
fl9|| . This fundamental feature holds for both uniform and non-uniform Bose systems, as 
can be seen, for example, from Eq. (6.38) of Hohenberg and Martin [[|. Thus by using 
the HFB single-particle matrix Green's function G\ to generate the density response func- 
tion by functional differentiation, one is effectively generating an "improved" single-particle 
spectrum which will be gapless. In particular, as discussed on p. 350-351 of Ref. [|J, the 
HFB Gi generates a density response function which is found to be identical to Beliaev's 
second-order single-particle spectrum ||. In the uniform case, this equivalence has been 
formally proven at all temperatures in a Bose-condensed gas by Cheung and Griffin (see 
Section IV of Ref. ||). The fact that this new spectrum is gapless follows from the fact that 
the Beliaev self-energies satisfy the HP relation in (^), with fi being given by the HFB 
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result in (g3J) . 

The above procedure shows that the excitations given by the Beliaev second-order ap- 
proximation is the correct generalization on the HFB excitation spectra. Such calculations 
involve evaluating the various bubble or polarization diagrams (see Eq. (2.30) of Ref. ||) 
which involve products of two HFB single-particle propagators G\. However, at temper- 
atures close enough to Tbe, it is adequate to evaluate these polarization bubble diagrams 
using a simple particle-like spectrum such as used in Refs. p|j20[]. 



Within a variational calculation, Bijlsma and Stoof have recently obtained (for a 
uniform system) the equivalent of the Popov approximation as defined above, but with g 
replaced by the many-body t-matrix calculated in the ladder diagram approximation. This 
extension and its relation to the full Beliaev second-order approximation will be discussed 



elsewhere 22 



IV. CONCLUSIONS 

Within the full Hartree-Fock-Bogoliubov approximation, we have given a simple deriva- 
tion of the coupled equations of motion for the condensate $(r) and its excited states, as 
given by (H) and (p3|) respectively. In contrast with previous work which was variational in 
nature |||3| , our derivation in Section II works directly with a mean-field approximation for 
the equation of motion for the condensate (<&) and non-condensate (ip) parts of the quantum 
field operator. This brings out the physics involved most clearly and has the advantage that 
one can make contact with the more general Green's function formulation [|]J| discussed in 
Section III. The self-consistent HFB has, of course, the advantage of being derivable vari- 
ationally and thus the total energy is minimized even if the HFB excitations have certain 
deficiencies. 

Because of the dynamical correlations induced by the Bose-broken symmetry, a consis- 
tent theory of excitations is surprisingly difficult to formulate even in a dilute Bose gas. 
This is already shown when one tries to improve the "gapless" Bogoliubov approximation 
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by working with the "conserving" Hartree-Fock-Bogoliubov (HFB) approximation which is 
not "gapless" . Hohenberg and Martin (HM) give a systematic discussion of this problem, 
introducing a classification of different approximations. In the present paper, we have em- 
phasized the usefulness of the HM classification when dealing with inhomogeneous atomic 
Bose gases in an external potential well |rj. We have used it to examine the single-particle 
spectrum predicted by the full HFB [g] as well as various simplifications of it which have 
been used in the literature PH|7| . In addition, the HFB has a special significance since 
it can be used to generate a density response function by functional differentiation which 
turns out to have the same spectrum as the Beliaev second-order approximation || for the 
single-particle Green's functions at all temperatures below Tbe- 

The reason that finding the "correct" excitation spectrum of a Bose-condensed gas is 
complicated is because it is difficult to simultaneously satisfy different requirements within 
a given approximation. The Hohenberg-Martin classification scheme we have used in Sec- 
tion III of this paper does not solve the problem of finding an approximation which is man- 
ifestly both gapless and consistent with conservation laws. An alternative approach which 
addresses this problem is based on the dielectric formalism (for a review, see Section 5.1 of 
Ref. [|n|)- This is a diagrammatic procedure which builds in the close connection between 



the single-particle and density-response functions so that their excitation spectra are iden- 
tical, with the conservation laws being incorporated through generalized Ward identities. 
This formalism |23] cuts across the HM classification and gives a procedure for choosing 



approximations which yield both a gapless elementary-excitation spectrum and response 
functions consistent with conservation laws (sum rules). In a future publication, we hope to 
discuss inhomogeneous Bose gases using this dielectric formalism approach ||23l , |T9| . 
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